Experimental and numerical signatures of dynamical crossover in orientationally 

disordered crystals 
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By means of NMR experiment and MD computer simulation we investigate the dynamical prop- 
erties of a chloroadamantane orientationally disordered crystal. We find a plastic-plastic dynamical 
transition at Tx — 330 K in the pico- nanosecond regime. It is interpreted as the rotational analogue 
of the Goldstein crossing temperature between quasi-free diffusion and activated regime predicted 
in liquids. Below T^, NMR experimental data are well described by a Frenkel model corresponding 
to a strongly anisotropic motion. At higher temperatures, a drastic deviation is observed toward 
quasi-isotropic rotational diffusion. Close to Tx, we observe that two-step relaxations emerge. An 
interpretation which is based on the present study of a specific heat anomaly detected by a recent 
calorimetric experiment is proposed. 



The elusive nature of the glass formation relates to 
the sharp rise of the transport coefhcicnts (viscosity or 
relaxation times) in a narrow temperature range above 
the calorimetric transition temperature Tg j^]. Some of 
the recent active investigations focus on the existence 
of a dynamical crossover which could occur above Tg in 
glass-forming liquids. The idealized version of the Mode 
CoupHng Theory(MCT) predicts that a dynamical de- 
coupling occurs at a critical temperature Tc > Tg It 
is however suspected that Tc corresponds to or is close 
to the temperature T^ at which the system also starts 
being sensitive to the energy landscape in its exploration 
of the configurational hyperspace ^-§1 . Motions in liq- 
uids would thus be dominate at low temperature by ac- 
tivated processes between potential energy minima while 
free diffusion occurs at high temperatures. Such dynam- 
ical crossover was already argued by Goldstein thirty 
years ago In the last few years a significant ef- 

fort has been undertaken to investigate experimentally 
and numerically the properties of glass-forming liquids 
in order to detect such dynamical changeover. However, 
most of the good glass-formers which allow easy inves- 
tigations of the supercooled metastable state such as o- 
terphenyl 1^ are molecular liquids where contributions 
due to translational, orientational and internal degrees 
of freedom (TDOF, ODOF) are inextricably mixed. It 
makes the properties of the energy landscape topography, 
MCT predictions or detection of dynamical transitions 
extremely hard to investigate. A clear visualization of the 
different processes at the molecular level of simple model 
compounds is needed. The discovery that some orienta- 
tionally disordered crystals exhibit the phenomenology of 
conventional glasses is therefore an exciting opportunity. 

Indeed, some molecular crystals present a plastic phase 
in which the molecular centers of mass show a perfect 
average crystalline translational order while the orienta- 
tions are dynamically disordered. Some plastic crystals, 
i.e glassy crystals can be undercooled and present 

many properties characteristic of conventional molecular- 



liquid glasses, both in the way the glass occurs (under- 
cooling) and the thermodynamic signatures. Glassy crys- 
tals which are the very rotational analogue of liquid glass- 
formers offer an excellent opportunity to focus on the role 
of the ODOF. In |jn|, Barrat et Al. investigated the ro- 
tational dynamics of a colliding hard needle model on a 
rigid lattice and have shown the emergence of a dynam- 
ical decoupling. Substituted adamantanes plastic crys- 
tals such as the cyanoadamantane |^2| noted CNa in the 
following are good experimental candidates and some of 
them exhibit a glass transition signature. In MD nu- 
merical simulations of a simple model inspired from the 
CNa molecular geometry allowed us to investigate both 
TDOF and ODOF properties. In this work a dynami- 
cal decoupling and two-step relaxations were found to be 
possible in a rotator phase. 

In the present letter, we address the possibility of 
finding such a behaviour in a real system. We have 
found chloroadamantane C10H15CI (Cla) to be a very 
favorable and relevant system for this investigation. It 
shows a plastic phase over a wide interval of tempera- 
ture [244 — 442] K (see Fig. |l|c) and a structure isomor- 
phous to CNa [Ol . The Cla molecule possesses a smaller 
substitute than does CNa. This gives rise to faster dy- 
namics in the plastic phase where rotational motions are 
also suspected to change in nature as it is reported from 
earlier incoherent quasielastic neutron scattering exper- 
iments Our results show for the first time, for a 
rotator phase, using experimental NMR and numerical 
MD simulations, the evidence of a dynamical crossover 
between quasi-free rotational diffusion and activated mo- 
tion occurring in the pico-nanosecond regime. It is also 
shown that this accident can explain the existence of a 
calorimetric anomaly (see Fig. ^) recently detected in 
Cla by Oguni et Al. 

NMR experiments were performed on Chloroadaman- 
tane on Bruker spectrometers ASXIOO and AMX 400. 

and ^'^C linewidths and relaxation times were mea- 
sured. The experiments were carried out on powder sam- 
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pies sealed in glass tubes for experiments or packed 
into zirconium rotors for ^^C MAS experiments. The 
spinning speed of the sample was 4 kHz. The sample 
temperature was regulated to within ±0.5 K by a Brucker 
BVT 2000 temperature controller. The spin lattice relax- 
ation times were measured using the inversion recovery 
pulse sequence, the recovery of the magnetization was 
always found to be exponential within experimental ac- 
curacy. The spin lattice Ti vs temperature is reported 
in Fig. |l|a. The main feature of this curve is the occur- 
rence of an abrupt change of the activation energy, from 
1400 K to 3400 K, at a temperature close to ~ 330 K 
which is far from both the melting temperature and the 
brittle to plastic transition. 
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FIG. 1. Experimental data as function of temperature in Cla 
plastic phase, a, the experimental proton relaxation time T\ and 
the dipolar order T^d- b, the ratio Ti (C4)/ri (C2) of the methy- 
lene relaxation times. The chloroadamantane molecule with the 
different carbons is displayed in the inset, c, Experimental heat ca- 
pacity where data have been extracted from Oguni et Al (Ref [ p^ ) . 

Slow diffusive motions of the molecules are generally 
observed in plastic phases. In the present study, their 
contribution to the spin-lattice relaxation time T\ was 
found to be negligible in the temperature range inves- 
tigated. This fact was particularly evidenced by mea- 
surements of the dipolar order Tid (see Fig. |l|a). The 
diffusional correlation time [jdiff — 30000 s at room 
temperature) was found to be strongly temperature de- 
pendent with a very high activation energy ~ 14000 K. 
Self diffusion has thus to be discarded as the origin of the 
break observed in Ti . Previous neutron studies have 
suggested that two types of motion occur in the rotator 
phase of Cla: tumbling of the molecular C3 axis among 
six fourfold crystallographic axes along the [100] direc- 



tions, and fast uniaxial rotation about the threefold axis. 
It was shown that the motion around the molecular axes 
is very fast (~ 10"^^ s at room temperature) with an ac- 
tivation energy of about 1240 K. Therefore, it takes place 
outside the time window where the dynamical crossover 
is detected at T^. Dipolar tumbling motions are thus 
mainly involved in the evolution of the spin-lattice re- 
laxation. Two separate temperature domains are clearly 
identified: below T^, both NMR experimental proton and 
carbon data are well described by a Frenkel model as- 
suming an Arrhenian behavior of tumbling and uniaxial 
rotations. Above T^, it was not possible to describe the 
evolution of Ti using the previous model. Hence, we ex- 
pect that some hypotheses of this model in particular the 
finite number of equilibrium positions become irrelevant 
at high temperatures. 

Anisotropy in the rotational motion can be precisely 
seen from the spin lattice relaxation times of the three 
different protonated carbons. Effectively the main relax- 
ation process of the methylene and methine carbons is 
the dipole-dipole interaction between the ^^C nuclei and 
their attached protons (see Fig. |l|b). Relaxation times 
for carbon were measured as a function of temperature 
in the rotator phase. The ratio Ti(C4)/Ti(C2) of the re- 
laxational times for the methylene (C2, C4) is displayed 
in Fig. A strongly anisotropic motion is found below 
Tx where the ratio of the carbon relaxational times is 
about 0.25. This latter value can be well reproduced by 
the previous Frenkel model. A drastic deviation from the 
anisotropic behaviour occurs at high temperatures where 
a ratio value of 0.7 is reached. Since free isotropic rota- 
tional diffusion processus would correspond to identical 
relaxational times for the methylene, this result reveals 
the influence of the residual crystalline field on the ro- 
tational dynamics. The experimental heat capacity ob- 
tained by Oguni et Al. |]1^, in the Cla plastic phase 
is reported in Fig. |l|c. The evolution of Cp shows an 
anomalous hump which was reported to resemble some- 
what to a Schottky anomaly. Owing to the nature of 
the Cla molecule which is mainly rigid, conformational 
changes have to be discarded for the origin of the Schot- 
tky anomaly. Clearly, the hump covers the entire tem- 
perature range where the dynamical crossover is seen in 
our NMR investigation. Its origin will be discussed in 
the light of the MD simulation results. 

MD simulations have been performed at 29 different 
temperatures from T = 220 to 500 K for a sample cor- 
responding to the crystalline fee rotator phase on a sim- 
ple model of Cla (see Fig. p|) A complete description 
of the model is given in ||l^,^8|. Owing to the very long 
MD runs, we succeeded in investigating truly equilibrated 
states of the system down to 220 K. The reorientational 
motions can be described by a set of correlation functions 
which are defined as: C;(t) = J2i=i,N i^i (^i W-A'i(O))) 
where P; is the Z-order Legendre polynomial and fli is the 
individual dipolar moment of the molecule i. The nature 
of dynamical changes and anisotropy can be checked by 
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computing Ci=i,2 functions which can also be related to 
the information obtained from NMR relaxation measure- 
ments. 

Fig. H shows the C2(t) time correlation function for 
all investigated temperatures above 220 K. Clearly, a 
two step relaxational behaviour classically seen in su- 
percooled molecular liquids is shown in our sim- 
ulations. When lowering the temperature, an intermedi- 
ate plateau region emerges which proves the presence of 
an orientational caging between neighboring molecules. 
This is the rotational analogue of the translational cage 
effect observed in liquids. This transient regime is fol- 
lowed by a slow a process which can be associated in our 
system, as checked for CNa , with large tumbling mo- 
tion between fourfold crystallographic directions where 
dipoles are preferably localized (see Fig. |^). A conse- 
quence is that the tumbling of one molecule is allowed 
only if an orientational rearrangement of its local neigh- 
bors occurs. This cooperative motion is clearly displayed 
using projections of the individual dipolar moments on 
one crystallographic plane in Fig. |^. 

(a) (b) 
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FIG. 2. a. Snapshots of the fl dipolar moments are projected on 
one xy plane at different instants in a run of 250 ps above (400 K) 
and below (220 K) Tx ~ 330 K. At T = 220 K, Some orientations 
are not populated since the system is not completely equilibrated 
over this duration, b, trajectory of the dipolar cosine angle 6 of a 
target molecule at the same temperatures. 

The Ti and T2 relaxational times are defined as the time 
it takes for their respective time correlation functions to 
decay e~^ of their initial values. They are displayed in 
Fig. ^a. Clearly, close to ~ 330 K, both n and T2 
relaxational times start diverging simultenously from the 
Arrhenian high temperature behaviour. Both high and 
low temperature apparent activation energies are found 



to be in good agreement with NMR experimental results. 
The nature of the rotational motions involved in this dy- 
namical change is given by the ratio ti/t2 in Fig. 
This latter evolves from a value equal to 3 correspond- 
ing to free small-step rotational diffusion to a value of 1 
classically associated to activated jump-like motion (see 
also Fig. H) |^,|2^. These results prove for the first time 
the existence of a plastic-plastic transition which can be 
interpreted as the rotational analogue of the Goldstein 
crossing temperature in liquids. It could also be associ- 
ated with a different process of exploring the energy land- 
scape topography above and below T^. As it has been 
recently proposed, dynamics of the system can be 

viewed as hopping processes between separated inherent 
structures at sufficiently low temperatures. 
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FIG. 3. Orientational correlation function C'2 as function of 
time t at different temperatures from 220 to 500 K. 

Specific heats have been directly extracted from the 
calculated enthalpy {Cp = dH/dT) and are shown in 
Fig. Clearly, a pronounced hump is seen which ex- 
tends over the low temperature part of the rotator phase 
below 350 K. It proves that our MD simulation captures 
the main features which are responsible to the anomaly 
found in calorimetric experiments by Oguni et Al. ]l6[ | . A 
good agreement with the experimental results is found for 
the area under the Cp hump and its temperature exten- 
sion while for the position of the maximum it is only fair. 
Obviously, the underlying absolute value of the specific 
heat is not reproduced in this computation owing to the 
simplicity of our model, where the internal vibrational 
contributions of the Cla molecule and the fast uniaxial 
rotations have not been taken into account. Thcrmody- 
namically, the anomalous hump requires an abrupt de- 
crease and a sharp change in curvature of the entropy vs 
temperature evolution. Fig. ^ shows that as the tem- 
perature is lowered the Cp hump begins to develop when 
dynamics change in nature. The strong localization of the 
dipole along the [100] directions (see Fig. ||) which devel- 
ops below Tx, is thus the prime suspect for the drastic 
change of the configurational entropy. 
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FIG. 4. a, orientational a-relaxation time ri (circle) and T2 
(square) obtained by simulation as function of the temperature T. 
The quantity plotted is T ln(Ti^2/''"o) where the time tq is obtained 
by fitting the relaxation times values for T = [400 — 500] K to 
an Arrhenius form. Deviation from the constant high-temperature 
value ~ 1418 K (dashed line) is observed around Tx ~ 330 K. Low 
temperature apparent activation energy is about 3109 K for ri and 
3325 K for T2 (solid line) in good agreement with NMR experimen- 
tal data, b, the ratio of the relaxational times ti/t2. c, calculated 
heat capacity as function of the temperature. 

In conclusion, convincing experimental and numeri- 
cal evidence for a dynamical crossover (T^) in the rota- 
tor phase of chloroadamantane have been obtained, very 
much in keeping with recent views on structural glass- 
formers. It is seen as a transition from quasi-free ro- 
tational diffusion of the molecular dipoles to activated 
geared tumblings at low temperatures. This rotational 
localization is shown numerically to give rise to a specific 
heat hump as was recently measured |l6| . The possible 
existence of such a calorimetric signature of a dynam- 
ical accident in liquids is highly debated and was un- 
successfully searched until now 2^ . Certainly, the 
underlying lattice and specificity of the orientational dis- 
order contribute to make dynamical changes in the sys- 
tem much sharper than in liquids. We now feel confi- 
dent that orientationally disordered crystals provide a 
valuable analog and alternative to the conventional liq- 
uid glass-formers. Results of the present study offer new 
interesting possibilities of testing the different theoretical 
approaches of the glass formation. 
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TERREG II program (Nord Pas de Calais/Kent). The 
authors wish to acknowledge the use of the facilities of 
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